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We report a detailed neutron scattering study of the ordering of spins and holes in oxygen-doped 
La2Ni04.i33- The single-crystal sample exhibits the same oxygen-interstitial order but better de- 
fined charge-stripe order than that studied previously in crystals with 5 — 0.125. In particu- 
lar, charge order is observed up to a temperature at least twice that of the magnetic transition, 
T m — 110.5 K. On cooling through T m , the wave vector e, equal to half the charge-stripe density 
within an MO2 layer, jumps discontinuously from | to 0.2944. It continues to decrease with further 
cooling, showing several lock-in transitions on the way down to low temperature. To explain the 
observed lock-ins, a model is proposed in which each charge stripe is centered on either a row of 
Ni or a row of O ions. The model is shown to be consistent with the Z-dependence of the magnetic 
peak intensities and with the relative intensities of the higher-order magnetic satellites. Analysis of 
the latter also provides evidence that the magnetic domain walls (charge stripes) are relatively nar- 
row. In combination with a recent study of magnetic-field-induced effects, we find that the charge 
stripes are all O-centered at T > T m , with a shift towards Ni centering at T < T m . Inferences 
concerning the competing interactions responsible for the the temperature dependence of e and the 
localization of charge within the stripes are discussed. 
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I. INTRODUCTION 

La2NiC>4 is a Mott-Hubbard insulator consisting of 
antifcrromagnetic MO2 planes alternating with La2C>2 
layers. The N1O2 planes can be doped with holes; 
however, contrary to conventional expectations, the 
material remains Jionmetallic up to quite large hole 
concentrations.BEm There is now considerable evidence 
that the insulating behavior occurs because the dopant- 
induced holes rfgwhto Pf 4s r themselves in periodically 
spaced stripesu'u'&Qffl'Ej'EJ These charge stripes run di- 
agonally relative to the square lattice defined by the 
Ni-O-Ni bonds. In the essentially undoped regions be- 
tween the stripes the Ni spins can order antiferromag- 
netically, w Vthijth l |' gjh ^ l | T *Sf f tripes acting as antiphase do- 
main walls .BEIj BHbIiI'EI This behavior represents an 
excellent example of "topological" doping in a correlated 
insulator £-3 

In comparing results obtained on samples doped both 
by Sr substitution and by addition of excess oxygen, 
several trends have become clear. In samples where 
both charge and spin orderings have been observed, the 
charge generally orders at a higher temperature thanj-the 



spins Eolld (The one apparent exception is the caseS of 
La2Ni04+,5 with S = 0.125. It will be shown in the 



present paper that, in this composition as well, the charge 
orders first.) Furthermore, both the charge and spin or- 
dering temperatures increase ag|BEaximately linearly with 
the net hole concentration n/j.QtJ'LS Also varying linearly 
with nh is the low-temperature spacing between charge 
stripes, corresponding to a hole density of roughly one 
per Ni site along a stripe. 

The observed low-temperature hole density within the 
stripes is consistent with-the mean-field calculations of 
Zaanen and Littlewood.cJ They evaluated an appropri- 
ate three-band Hubbard model at zero temperature and 
found that electron-phonon couplings tend to reinforce 
the stability of charged domain walls in an antiferromag- 
netic background. The numerical calculations yielded a 
relatively narrow charge stripe centered on a diagonal row 
of Ni atoms. Zaanen and Littlewood pointed out that 
the domain wall states form a half-filled one-dimensional 
(ID) band, which is unstable to a Peierls distortion. Such 
a distortion would create a gap at the Fermi level, thus 
explaining the insulating behavior. 

Despite the apparent success of the Hartree-Fock 
calculationsEj in describing the low-temperature stripe 
structure, such an approach runs into difficulties at high 
temperature. In the absence of stripe order, a Hartree- 
Fock analysis of the doped nickelate w^uk^j^d-Sr^jnetal- 
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lie state, whereas numerous studies 



have 
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shown that even above the charge-ordering temperature 
the optical conductivity is dominated by midinfrared ab- 
sorption bands, with no ui = Drude component. Fur- 
thermore, neutron diffraction studies of the ordered state 
have shown that iruppme cases the stripe spacing changes 
with temperature,ErEI implying a variation in the charge 
density within a stripe. Thus, one cannot rely on a 
unique stripe charge density to explain the insulating be- 
havior. 

It is clear that there are mysteries remaining to be 
unraveled. For example, what is the nature of the 
temperature-dependent variation of the stripe spacing? 
How narrow are the charge stripes, and to what extent 
are they pinned to the lattice? Are charge stripes ac- 
tually centered on Ni rows, or are other configurations 
possible? To address these issues, a well-ordered sample 
is required. 

Much of the work cited above has focussed on Sr- 
doped La2NiC>4. That system has the advantage of al- 
lowing a continuous variation of the hole concentration 
by dopants whose random positions are fixed at rela- 
tively high temperature. The down side is that the dis- 
order of the dopants appear^a -prevent the development 
of long-range stripe order .QLjtrEJ The alternative is to 
study oxygen-doped La2NiC>4 +( 5. The constraint here is 
that one must work with the phases that nature pro- 
vides: ordered phases occur only at special values of 6, 
and for an.- arbit. i jaipf ii vaVue <? f S one may observe phase 
separation.Eil3c3MESE3'E§EZl On the positive side, long- 
range order of intcrstitials and charge stripes coincidft4JtA 
composition originally characterized as 5 = 0.125.li&tll'Ej 

Analysis of neutron diffraction measurements on one 
particular 5 = 0.125 crystal! indicated that the ordered 
interstitial structure actually corresponds to an ideal oxy- 
gen excess of 6 = jg = 0.133. Thus, the 6 = 0.125 com- 
position is oxygen deficient relative to the ideal phase. 
As will be explained below, we have prepared a new crys- 
tal that appears to correspond much more closely to the 
ideal value of 8. The magnetic Bragg peaks measured 
on this crystal arepSharper than those observed with the 
5 = 0.125 crystal,™! thus allowing better sensitivity to 
the features of interest. We have recently reported results 
on magnetic-field-induced staggered magnetizatioi£j and 
on spin excitations in the stripe-ordered phaseEJ mea- 
sured on the same crystal. In the present paper we con- 
centrate on various aspects of the static stripe order in 
zero field. . 

In earlier studiescrEa of S = 0.125 crystals, an appar- 
ently discontinuous jump from zero to finite intensity of 
both the magnetic and charge-order superlattice peaks 
was observed on cooling below T m = 110.5 K. With the 
S = 0. 133 crystal we observe a similar discontinuous jump 
in the intensities, but we find that charge order clearly 
survives at T > T m with a jump in the stripe spacing 
to a value commensurate with a second harmonic of the 
interstitial order. In the latter regime the charge-order 
peak intensity decreases exponentially with temperature, 
so that no clear transition to the disordered state has 



been observed. 

To obtain a more precise measure of the incommen- 
surate splitting e of the magnetic peaks at T < T m , 
we measured the temperature dependence of the posi- 
tion and intensity of a third harmonic peak. Where pre- 
vious workQ had provided weak evidence that e, which 
varies with temperature, tends to show lock-in behav- 
ior at certain rational fractions, the present results show 
clear plateaus at e = § (98 K< T <102 K) and e = | 
(75 K< T <92 K). The lock-ins demonstrate that the 
charge stripes are coupled to the lattice. To explain 
the continuous distribution of e values found at other 
temperatures, an analogy is made with the problem of 
finding the ground state arrangement of integral charges 
on a one-di m e ns ional lattice at arbitrary filling. From 
the solutiontiirc3 to the latter problem, it appears that 
a devil's staircase of ordered phases should be possible; 
however, due to entropy, long-range order is observed 
only for phases with relatively short periods. The fact 
that e varies with temperature indicates that competing 
interactions must be involved in determining the order .□ 

Coupling to the lattice is also evident in the way that 
the stripes align themselves from one layer to the next. 
If there were no interaction with the lattice (as was as- 
sumed previouslyu) , then the stripes within a layer would 
be uniformly spaced, and their positions would be stag- 
gered from one layer to the next in order to minimize the 
Coulomb energy. Such a symmetric stacking of the layers 
of stripes can lead to forbidden superlattice peaks corre- 
sponding to the charge order, similar to the situation 
for nuclear reflections due to the average body-centered 
stacking of the NiC>2 layers. However, the pinning of the 
charge stripes to the lattice means that the shift of the 
stripe pattern from one layer to the next can only occur in 
increments of the lattice spacing. As a consequence, the 
intensities of the magnetic peaks oscillate as a function 
of the component of the momentum transfer that is per- 
pendicular to the planes. We show that the modulation 
of the intensities can be modeled quantitatively in terms 
of the stacking structure factor when a small amount of 
stacking disorder is taken into account. It follows that 
the missing charge-order peaks are not forbidden, but 
just extremely weak. 

To characterize the widths of the charge stripes, we 
have measured the higher harmonic magnetic Bragg 
peaks. Narrow domain walls give squared-off magnetic 
domains, and deviations from a sinusoidal modulation 
are reflected in finite Fourier components beyond the first 
harmonic. Another contribution to the higher harmonics 
can come from an alternation of the stripe spacing as- 
sociated with the pinning to the lattice. We show that 
although the intensities of the higher harmonics are all 
quite weak, with the strongest (the third) being < 2% of 
the first, they are consistent with narrow charge stripes. 
Furthermore, the weakness of the harmonics in the case 
of e = | indicates that stripes can be centered both on 
rows of Ni atoms and on rows of O atoms. Assuming that 
the stripes are Ni-centered at low temperature, the tem- 
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perature dependence of e implies that the stripes become 
O-centered at high temperature. This leads to a pre- 
diction of ferrimagnetic correlations in the stripe-ordered 
regime at T > T m , which was recently confirmed by mea- 
surements in a magnetic fieldE3 

The rest of the paper is organized as follows. The 
next section contains a description of experimental pro- 
cedures, and it is followed by a brief review of nomen- 
clature and notation. The temperature dependences of 
the intensities and positions of superlattice peaks are pre- 
sented and analyzed in Sec. IV. The results relevant to 
interlayer stacking of stripes are treated in Sec. V. In 
Sec. VI the intensities of the magnetic higher harmonic 
satellites are analyzed. The significance of these results 
is discussed in Sec. VII, and a brief summary is given in 
Sec. VIII. 



II. EXPERIMENTAL PROCEDURES 

The single-crystal sample used in the present study 
was grown-by radio frequency skull melting as described 
elsewhere.cJ The oxygen concentration was selected by 
annealing at 464° C for 5 days in flowing O2 (1 bar) 
followed by a quench to room temperature. The anneal- 
ing conditions were chosen to give a nominal-^ of 0.150, 
based on earlier work on the phase diagram.c-3 An initial 
neutron diffraction study of this crystal in 1993 indicated 
that the dominant oxygen-ordered phase is the same as 
that found previouslyO in a crystal with S = 0.125. There 
were also weak diffraction peaks indicative of a second 
phase with a different <5; however, these were absent ~ 1.5 
years later when the present study was initiated. (One 
fact that might be relevant to the disappearance of the 
second phase is that the crystal was stored in a He at- 
mosphere, which is relatively reducing, rather than in 
air.) The room temperature lattice parameters obtained 
at that later time are a = 5.461 A and c = 12.674 A. 
Making use of the calibration established by the earlier 
X-ray diffraction study,Eil the c/a ratio of 2.3208 indi- 
cates 5 = 0.136± 0.005, consistent with the ideal value of 

-£r for the model ordered interstitial structure identified^ 
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in the case of S = 0.125. 

The neutron diffraction measurements were performed 
on the H4M, H8 and H9A triple-axis spectrometers at 
the High Flux Beam Reactor located at Brookhaven Na- 
tional Laboratory. At H9A, which is on the cold source, 
neutrons of energy 5.0 meV were selected. Effective hori- 
zontal collimations of 40'-40'-20'-40' (from reactor to de- 
tector) were used together with a cold Be filter to el- 
liminate neutrons at shorter harmonic wavelengths. At 
H4M and H8 a neutron energy of 14.7 meV, horizon- 
tal collimations of 40'-40'-40'-40' or 40'-40'-80'-80' and a 
pyrolytic graphite (PG) filter were used. At all spectrom- 
eters the neutrons were monochromatized and analyzed 
by PG crystals set for the (002) reflection. 

The cold neutrons at H9A were used to study the 



weak higher-order satellites of the magnetic scattering 
peaks, which are strongest at small momentum transfer 
Q. The excellent signal to background ratio of H9A to- 
gether with the larger scattering angles in the case of cold 
neutrons compared to higher energy neutrons enabled the 
determination of third-order satellite positions with high 
precision. The higher energy neutrons at H4M and H8 
were required to measure the nuclear scattering associ- 
ated with lattice displacements due to charge order. Ex- 
cept for the third-harmonic measurements, where scans 
were taken along Q = (h, 0, 1), integrated intensities were 
determined from 6 — 29 scans. 

All measurements were performed with the crystal ori- 
ented so that the (h, 0, 1) zone of reciprocal space was in 
the horizontal scattering plane. (The notation and choice 
of unit cell are explained in the following section.) The 
aligned crystal was mounted in an Al can filled with He 
exchange gas that was attached to the cold finger of a 
Displex closed-cycle refrigerator. The temperature was 
measured with a Si diode and a Pt resistance thermome- 
ter to an accuracy of better than ±1 K. 

The ordering kinetics of the oxygen interstitials can 
be rather slow, and the cooling rate through the or- 
dering transition can affect whether long-range order is 
achieved. For example, extremely slaffi cooling is required 
to optimize the order in the casd~3 of S = ClOS^r-apd 
rapid cooling also had measurable effects in studiesH'tL3 of 
5 = 0.125. In the present case, no significant effects of 
cooling rate on interstitial order were observed. On the 
other hand, the cooling and heating history for T < T m 
has an influence on the relative fractions of coexisting 
magnetic phases, the width of the temperature regions 
where the magnetic wave vector is locked-in to rational 
fractions, and the detailed spectral weight of higher order 
magnetic satellites and their peak widths. Therefore, the 
temperature dependence was studied during slow, careful 
warming after an initial cool down to 10 K. 

III. NOTATION 

In order to describe and explain the data it is first nec- 
essary to review the structureEl and explain our notation. 
Within the resolution of our instrument, the fundamen- 
tal reflections of this crystal indicate a tetragonal (HTT) 
symmetry (space group IA/mmm) down to at least 10 K, 
whereas the structural superlattice peaks suggest an or- 
thorhombic symmetry. To describe the diffraction peaks 
we will use an indexing based on an orthorhombic cell 
(or pseudo-tetragonal F4/mmm) with a unit cell size 
\J~2at x V2at x ct relative to the basic body-centered 
tetragonal cell. Note that in the body-centered tetrago- 
nal cell the a and b axes are parallel to the Ni-0 bonds 
within the MO2 planes, whereas in the larger cell that 
we use the axes are rotated from the bonding directions 
by 45°. 

Four different types of superlattice peaks can be found 
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in this sample. Two of them correspond to the ordering 
of the oxygen interstitials and the other two involve or- 
dering of the doped holes and the Ni spins. Each set of 
peaks is characterized by a wave vector gj such that Qsl, 
the position of a superlattice peak of type j, is given by 



g„ £ = (ne,0,0). 



(6) 



Q S L = G± Sj , 



(1) 



where G is a reciprocal lattice vector corresponding to 
the average pseudo-tetragonal unit cell. The oxygen or- 
der is characterized by two wave vectors, 



goi 



1 



and 



g02 = 



,0,1 



4 4 
5'5 



(2) 



(3) 



where wave vector components are specified in reciprocal 
lattice units of (2w/a, 27r/a, 2w/c). The twinning associ- 
ated with the interstitial order is described in Ref. |^. 

To describe the superlattice peaks associated with 
stripe order it is necessary to expand the set of recipro- 
cal lattice vectors considered in Eq. ([]]). The symmetry 
relationship between the two NiOa planes per F4/mmm 
unit cell leads to forbidden reflections. For example, for 
reciprocal lattice vectors of the type (H, 0, L) one has the 
requirements that both H and L be even integers. We 
will show in Sec. V that the stripe patterns in neighbor- 
ing layers are not related by symmetry, and as a result 
the relevant reciprocal lattice vectors include the case L 
odd. 

The positions of the magnetic and charge-order super- 
lattice peaks within the (h, 0, 1) zone are illustrated in 
Fig. 1(a). The charge order is determined by 



g2e 



(26,0,0), 



(4) 



with the average distance between domain walls in real 
space equal to a/2e. The magnetic order involves two 
modulation wave vectors. Simple antifcrromagnetic or- 
der would be described by g = Qaf = (1,0,0). Stripe 
order results in a second modulation 



(e,0,0), 



(5) 



so that the net wave vector is g = QAF±ge- In real space 
the spin structure consists of locally antiferromagnetic 
order with an overall modulation period a/e, twice that of 
the charge period, reflecting the fact that the hole stripes 
act as antiphase boundaries for the spin structure. The 
idealized real-space spin and charge structure for ( = \ 
is presented in Fig. |](b). 

Of course, g e and g 2e represent only the dominant si- 
nusoidal components of the magnetic and charge-order 
modulations, respectively. Deviations from sinusoidal 
modulations will result in superlattice peaks correspond- 
ing to higher Fourier components of the order. The mod- 
ulation wave vector for the n th harmonic is simply 



For the magnetic scattering, only harmonics with n odd 
will appear, split about Qaf- (The absence of magnetic 
harmonics with n even follows directly from the antiphase 
nature of the order.) Harmonics with n even correspond 
to charge order. 



IV. TEMPERATURE DEPENDENCE OF STRIPE 
ORDER 

A. Experimental Results 

The overall temperature dependence of the stripe or- 
der is best characterized by looking at first and second 
harmonic peaks. Figure presents results for a mag- 
netic superlattice peak at (3 + e, 0, 1) and a charge-order 
peak at (4 — 2e,0, 1). The upper panel shows the inte- 
grated peak intensity normalized with respect to 10 K, 
while the lower gives the peak position h measured in 
a [h, 0,0] scan. As found previouslyJ§E3 cooling through 
T m = 110.5 K results in a discontinuous jump in inten- 
sity of both the magnetic and charge-order peaks, with 
a concomitant temperature dependence of the peak posi- 
tions. The new feature that we have to report is a peak at 
e = 3 for T > T m . Scattering at this position could come 
from either magnetic or charge order (or both); however, 
the Q dependence of the structure factor indicates that it 
corresponds uniquely to charge order. (The intensity of 
( = \ peaks disappears at small Q, where the magnetic 
form factor is largest.) 

The discontinuity of the transition at T m is illustrated 
in Fig. At a temperature 0.5 K below the transition 
we observe one peak due to magnetic order and another 
corresponding to charge order. Right at T m a third peak 
appears at e = i, coexisting with the previous two. At 
0.5 K above T m only the e = | charge-order peak sur- 
vives. In contrast to the behavior at T < T m , the position 
of the peak remains locked at e = ^ while its intensity 
gradually decreases with increasing temperature. 



The temperature dependence of the e = | peak is pre- 
sented on a logarithmic scale in the inset of Fig. |^. The 
linear variation of the intensity for T > 150 K is analo- 
gous to a Debye- Waller like behavior with 



j^ e -2T/T 



(7) 



and Tq = 67 ± 5 K. Such a temperature dependence is in 
sharp contrast to the usual critical behavior associated 
with an order-disorder transition. In connection with 
this, it is useful to note that 



(8) 



? 2£ = ( -,0,0] =2goi. 



so that there is a type of commensurability between the 
charge stripes in the planes and the interstitial order 
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along one direction of the lattice. The Debye- Waller-like 
decay of the intensity suggests that the charge correla- 
tions are fluctuating about an average ordered configura- 
tion determined by the ordered interstitial oxygens. It is 
possible that the charge order does not entirely disappear 
until the interstitials disorder. 

The ordering of the spins at T m is associated with 
jumps both in the charge-order intensity and in e. 
PreviouslyO it was found that the superlattice peak inten- 
sities show simple power-law correlations with the quan- 
tity 

q =l~e. (9) 

From the measurements between T m and 80 K, we find 
for the magnetic-peak intensity that 

imag ~ aq 2 , (10) 

as illustrated in Fig. || This result is completely empir- 
ical, as there is no theoretical prediction for the connec- 
tion between the intensity and wavevector. 

A prediction does exist for the relationship between the 
intensities of the magnetic and charge order peaks below 
T m . Zachar, Kivelson, and Emeryo have solved an ap- 
propriate Landau model for coupled charge- and spin- 
density-wave order parameters, p and S, respectively. 
They find that when the charge is the first to order, it is 
possible to have a first-order transition when the spins or- 
der, as we observe. At temperatures below the magnetic 
transition, the relationship between the order parameters 
is predicted to be, to lowest order, 

P -po = A\S\ 2 , (11) 

where A is a constant and po denotes the magnitude of 
p just above the magnetic transition temperature. Since 
the intensities of the magnetic and charge-order peaks 
are related to the order parameters by 7 ma g ~ |>5| 2 and 
J ch ~ \p\ 2 , Eq. (jll|) implies that, for T < T m , 

[7 ch (T)] 1/2 - [l ch (T+)] 1/2 = A'/ mag (T) (12) 



where I c h(T+) is the intensity of the charge-order peak at 
a temperature just above T m . Keeping in mind that the 
model does not describe the temperature dependence of 
the wave vector, so that there may be missing correction 
terms, we show in Fig. || that this relationship is at least 
plausible near the transition. 

The step-like structures in the lower panel of Fig. || 
are suggestive of lock-in transitions. To obtain more pre- 
cise measurements of e versus temperature, we studied 
a third-harmonic magnetic peak, (1 — 3e,0, 1). This re- 
flection occurs at a relatively low scattering angle, and 
a scan through the peak along the [h, 0, 0] direction is 
nearly transverse to Q. For such conditions the reso- 
lution in h is excellent. In our case the resolution was 



limited by the sample mosaic (« 0.7°). At many tem- 
peratures the peak shape and width are limited by the 
asymmetric mosaic distribution of the sample. Represen- 
tative scans are shown in Fig. ^. At other temperatures 
we observed either a single broadened peak, or the coex- 
istence of a broadened peak and a mosaic-limited one. 

A summary of the results is given in Fig. @. In all pan- 
els, the dominant component is represented by the filled 
symbols and the secondary component by the open ones. 
(The exception to this rule occurs in the range between 
50 K and 70 K, where the integrated intensity of the 
broader second component is actually greater than that 
of the sharper peak.) The top panel shows the integrated 
intensity, with the total intensity from both components 
(dashed line) normalized to one at 10 K. The middle 
panel and inset show e. The bottom panel gives the peak 
width, divided by 3 (in units of 2n/a), so that it can be 
compared directly with the middle panel. 

The temperature dependence of e is rather interesting. 
Although it was measured on warming, it is more conve- 
nient to discuss the behavior starting from T m and con- 
sidering the changes as the temperature decreases. (Hys- 
teresis effects should be small as long as one cools suffi- 
ciently slowly.) Initially, e jumps from | to 0.2944 w ^. 
From there it decreases linearly with temperature until 
it hits a plateau at 0.286 = |. (There is also a hint of 
a plateau at 0.293 = although the stability of this 
particular value is surprising, as discussed below.) The | 
plateau extends from 102 K down to 98 K. Below this, e 
jumps to another region of linear decrease, until it locks 
into a value of 0.2778 = jg near 93 K. From this point 
down to 10 K a ^ component is always present; however, 
the entire system is locked into this value only down to 
~ 75 K. Near 70 K the third harmonic splits into two 
peaks: a sharp one at and a broad one centered near 
0.275. The integrated intensity of the broad peak domi- 
nates down to ~ 50 K, below which most of the intensity 
shifts back, into the sharp jg peak. In contrast to what 
was foundQ with the 6 = 0.125 crystal, we have not ob- 
served any component centered at jy = 0.2727. 

B. Analysis 

In order to make sensep-of e(T), we can start with 
the Zaanen and Littlewoodll3 prediction of charge stripes 
with a hole density of one per Ni site. Assuming that 
all dopant-induced holes go into O 2p states lying within 
the MO2 layers, one would then expect a domain wall 
density of e = 26 per row of nickel atoms. In our case 
this gives e — ^ — 0.267. It appears that the measured 
e is approaching this value as the temperature decreases, 
but it never quite gets there. If we assume that all of the 
holes are indeed ordered in the stripes, then the hole den- 
sity (per Ni site) within a charge stripe, n s , is given by 
n s = 2S/e. Even if all holes remain confined to the stripes 
as the temperature is raised, the fact that e increases im- 
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plies that n s is decreasing. At 80 K we have n s = 0.960, 
at 100 K n s = 0.933, and at T > T m n s = 0.800. 

The existence of plateaus in e(T) indicates that the 
stripe order can couple to the lattice; however, the al- 
most continuous variation of e between plateaus shows 
that the stripe density is not restricted to a few special 
values. To describe the possible stripe arrangements, we 
need consider only the one-dimensional ordering along 
the charge-modulation direction [see Fig. 0(b)]. Within 
this model, let each row of Ni or O atoms be represented 
by a single site. Thus, there are two Ni sites per unit cell, 
separated by ^a, with an O between each pair of nick- 
els. The stripes may be pinned to particular sites in the 
lattice, either by electron-phonon coupling or Coulomb 
interaction with the ordered oxygen interstitials. The in- 
formation on the particular pinning sites is contained in 
the phase of the scattered neutron waves, which cannot 
be detected in our single-scattering measurements. To 
be definite, we will assume that a stripe can be pinned 
to either a Ni or an O row. (It will be shown below that 
this assumption leads to testable predictions.) Although 
it is not necessary for the present discussion, we will also 
assume that the charge stripes are quite narrow. (The 
actual charge distribution can be probed by analyzing 
the intensities of the higher-harmonic superlattice peaks, 
as will be discussed in Sec. VI.) For a Ni-centered stripe, 
we assume that there is no static moment on the sita 
where the hole is centered, consistent with experiment ,El 
but neighboring nickels have full spins. For an O-centered 
stripe, we assume that there is no static moment associ- 
ated with the hole spin. 

Since the stripes are charged, they will repel each 
other. As a result, the stripes will arrange themselves 
so as to maintain the maximum possible spacing, with 
the constraint that each stripe is centered on a Ni or O 
site. To understand the possible configurations, we can 
consider the equivalent problem of finding the ground 
state arrangement of integral charges on a one dimen- 
sional latticer-at arbitrary filling, which has been solved 
by Hubbard, PH hy Pokrovsky and Uimin,E3 and by Bak 
and Bruinsma.E3 For any given filling, there will be at 
most two stripe spacings, which differ by the minimum 
increment in stripe position. In our case, this increment 
is jo, corresponding to a shift from a Ni to an O site 
(or vice versa). The energy is minimized when the two 
spacing units are ordered in a periodic fashion. The most 
stable configurations are those with a single spacing. In 
our region of interest, the single-spacing phases occur at 
e = j, |, and ^, with stripe separations of 2a, ja, and 
|a, respectively. 

We can represent these different phases using a nota- 
tion in which an up or down arrow represents a Ni spin, 
a dot represents an O site, and an open circle represents 
a stripe. Thus, the configuration corresponding to e = \ 
with Ni-centered stripes looks like 




(13) 



Note that each spin domain has a net moment, but that 
the moments of neighboring domains are antiparallel. 
Shifting the stripes onto O sites gives 



4 it 



(14) 



where now we can associate moments with the domain 
walls. 

For e = |, we get equal numbers of Ni-centered and 
O-centered stripes: 

o-T-4-T°J-4-T-°-4-T-W-T-4-° (15) 

If we drop the requirement of a single spacing, we can 
have either purely Ni-centered stripes: 

°-T_4_T-°-4-T-°-4_T_4-°-T-4-° (16) 

or purely O-centered stripes: 

W-T-4-W-4-W-4-T-W-T-4 (17) 

ii- t 1r 4 

In the single-spacing phase e = i, Ni-centered stripes 
give 



°-W-°-M-o 
while O-centered stripes yield 




(18) 



(19) 



For our S = 0.133 crystal, the magnetic order vanishes 
when e shifts to ^; nevertheless, inelastic scattering mea- 
surements show that the spins remain correlated.E3 In the 
configuration with O-centered stripes, the domain walls 
have moments that are in phase with each other. As will 
be discussed further in Sec. VI, one might expect the mo- 
ments of Ni ions immediately adjacent to an O-centered 
domain wall to be somewhat reduced compared to those 
in the middle of a magnetic domain. If the modulation of 
the moments defining the antiphase domains were purely 
sinusoidal, then there would be no net moment in the 
system; however, any deviation from a pure sine wave 
would lead to fcrrimagnetism. In contrast, the spins all 
form antiparallel pairs when the stripes are centered on 
Ni rows. Recently we demonstrated^ that one can in- 
duce staggered magnetic order in the e = | phase by 
applying a uniform magnetic field. This result clearly 



G 



indicates that the stripes in the high-temperature phase 
are centered on O rows. 

So far we have considered only phases with uniform 
stripe spacing, whereas the real system appears to pass 
through a range of more complicated phases. While we 
have emphasized the possibility of two types of charge 
stripes, it is important to note that the observed values 
of e can be explained with a single type of stripe. In that 
case, the relevant phases of uniform spacing are limited 
to e = j and -|. To describe the possible intervening 
phases it is convenient to introduce a simple notation. 
First note that for e = j, there is 1 spin period in 4 lattice 
spacings, or, equivalently, 1 charge stripe per 4 rows of Ni 
atoms (with 2 rows per fundamental unit cell). It follows 
that this phase can be denoted in terms of the real-space 
repeat unit of the charge stripes as (4) . Similarly, e = | 
corresponds to (3). The configuration (43) represents 
alternating stripe spacings of 4 Ni rows and 3 Ni rows. 
This corresponds to 2 stripes per 7 Ni rows, or e = =. An 
arbitrary configuration will consist of n 4-row spacings 
and m 3-row spacings, giving)!] 



ri + m 
An + 3m 



(20) 



For such an arbitrary configuration, the ground-state en- 
ergy is minimized by arranging the .stripes in a pattern 
that is as close to uniform as possible.Eill2lSpme examples 
are shown in Table |. Bak and BruinsmaEj have demon- 
strated that the possible ground states form a complete 
Devil's staircase. 

Of course, if two types of stripes are considered, then 
we must include a possible stripe spacing of | Ni rows. 
Note that to be commensurate with the Ni-row spacing 
such units must occur in pairs. Let the number of pairs 
of ^-spacings per real space repeat unit be p. Then, in 
the range 4 < e < 1, the values of e are given by 



n + 2p 
4n + 7p ' 



(21) 



while, in the range j < e < ^, the appropriate formula 
is 

< 22 » 

Examples are shown in Table |l[ 

One can distinguish between models with either one 
or two types of stripes by analyzing the intensities of 
higher harmonic satellites. Such an analysis is discussed 
in Sec. VI. The proposed pinning of the stripes to par- 
ticular lattice rows has implications for the layer to layer 
stacking of the stripes, as we discuss in the next section. 

V. STACKING OF STRIPES 

To minimize the Coulomb repulsion between charge 
stripes in neighboring planes, the stripes should run par- 
allel to each other with a staggered alignment. Such 



a body-centered type arrangement would maximize the 
spacing between nearest neighbors; however, pinning of 
the stripes to the lattice would prevent such an ideal stag- 
gering. Within our model, the stripe pattern should only 
be able to shift by integral (or possibly half integral) units 
of the lattice spacing a. We can distinguish between these 
two possibilities, symmetric versus integral-shift stack- 
ing, by analyzing the I dependence of the superlattice 
intensities. 

Consider first the case of a symmetric, body-centered 
stacking of the stripes. As discussed in the previous sec- 
tion, a given value of e can be written as a ratio of two 
integers, r/s. The commensurate period of the magnetic 
structure is then sa, and the period of the charge struc- 
ture is sa/2. The displacement of the stripe pattern from 
one layer to the next is then (|;S,0, ^), in lattice units. 
The structure factor describing such a stacking is 



F = 1 



(23) 



where the intensity is proportional to \F\ 2 . For a charge- 
order reflection of the type (2n + 2e, 0, 1), where n is an 
integer, \F\ 2 is either 4 or 0. In the case of e = j, one 
finds that |_F| 2 is equal to for I equal to an even inte- 
ger and 4 for odd I, independent of n. For the general 
case of e = r/s, however, the value of |F| 2 depends not 
only on I, but also on r, s, and n. There is a similar 
variability for magnetic reflections. To illustrate this, we 
have listed values of |F| 2 for a number of different cases 
in Table [II. Note in particular that the body-centered- 
stacking model implies zero intensity for the (4 — 2e, 0, 1) 
reflection when e = |, which is clearly inconsistent with 
the results shown in Fig. |[ (The tabulated values of \F\ 2 
for e = i have implications for the interpretation of mag- 
netic .and charge-order scattering as reported in a recent 
studyB of Lai. 67 Sr .33NiO 4 .) 

Alternatively, if the stripes are pinned to the lattice, 
then the stripe order in neighboring layers is described 
by the basis vectors (0,0,0) and (ma,0, |), where m is 
an integer number. Taking the simplest case of m = 1, 
one finds that, for the n th harmonic superlattice peak, 



|F n | 2 = 2 [1 + (-l) z cos(27me)] 



(24) 



This formula indicates that, for fixed h, the intensities 
of magnetic peaks should oscillate as a function of Z. It 
turns out that this is just what is observed. 

We measured integrated intensities for 10 first- 
harmonic peaks at (1 ± e, 0,1), with I = to 4, at the 
80 K and 100 K plateaus, corresponding to £ = ^ and 
b, respectively. To describe the experimental results, it is 
helpful to express the integrated intensities of the n = 1 
magnetic reflections as 



h(l) 



AflF^ 2 , 



(25) 



where A is a normalization constant, and we use the 
parametrization / = e _< ^ / 2<T to approximate the Ni 
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magnetic form factor as well as the Debye- Waller fac- 
tor. The parameter a was determined by a least-squares 
fit to the data. The values of |-Fi| 2 extracted from the 
80 K measurements are shown in Fig. The dominant 
I dependence of these results is described by the formula 



\F n \ 2 =A'[l + (-l) l B n ], 



(26) 



single layer), we can calculate the relative intensities of 
even-l peaks compared to odd-Z peaks. Using our results 
for e = A gives a ratio of 2%, while for e = I we obtain 
8%. These relative intensities are still quite difficult to 
detect. For e = i we calculate a ratio of 33%, which 
would be detectable except that the peak position over- 
laps with the second harmonic of goi- 



where we have kept the notation general by labelling F 
and B with n, and the scale factor A' is meant to remind 
us that we have not determined from the experimental 
intensities the absolute magnitude of |Fi| 2 . The values 
of B characterizing the measurements at 80 K and 100 K 
are listed in Table [V . Equation predicts that B\ = 
cos(27re), which yields B x = -0.174 at 80 K (e = 
and B x = -0.223 at 100 K (e = f ). 

The model allowing an interlayer offset of the stripe 
pattern by one lattice spacing gives qualitative agreement 
with experiment, although the calculated magnitude of 
the parameter B\ appears to be systematically too large 
by 30%. We can get a better fit to the data if we allow 
for some disorder in the stacking in terms of different 
offsets. The Coulomb energy between charge stripes in 
neighboring layers is not the same for all possible inte- 
gral offsets. If we allow for two possible offsets, a and 
ma, with relative probabilities I — a and a, respectively, 
then, assuming a random distribution of the offsets, the 
structure factor becomes 



F n = l+ [(1 



-in I 



(27) 



In the yg phase, the Coulomb energy is minimized for 
m = 1, and also for m = 8. We fit the measured intensi- 
ties with this model, varying the parameters a and a to 
obtain the results shown in Fig. ||. The parameter values 
and the weighted reliability factor R w are listed in Ta- 
ble [v| For | , the relevant offsets are m — 1 and 3 (with 
a shift by 2.5a being equivalent in energy to m = 1). 
Again, the fitting results are listed in Table 0. The cal- 
culated values of B\ for both cases are given in Table 0. 

One can see from Eqs. ( pi[ ) and (27) that the relative 
size of the Z-dependent intensity oscillation varies with n, 
the harmonic order. Thus, another test of our model is 
to apply it to another value of n. We have parametrized 
measurements of (1 — 3e, 0, 1) peaks (I = 1 to 4) in terms 
of the coefficient B3, and the results for T = 80 K and 
100 K are listed in Table [V. We then calculated values 



of B 3 using the parameter values already determined in 
the fits to the n = 1 intensities. The calculated values 
are in very good agreement with the measurements. For 
reference, we have also included in Table IV calculated 
values of B2 and B4. 

Before leaving this section, it is of interest to consider 
the consequences of our model for the n — 2 charge-order 
peaks. Measurements have not detected any intensity 
at positions (2n ± 2e, 0, 1) with I even. Considering just 
the stacking contribution (that is, ignoring the structure 
factor due to the atomic displacement pattern within a 



VI. MAGNETIC HARMONIC SATELLITE 
SPECTRUM 



Information on deviations from a sinusoidal modula- 
tion is contained in the relative intensities of the higher 
harmonic superlattice peaks. In particular, we are inter- 
ested in the degree to which the spacing of the stripes 
deviates from perfect periodicity as well as the degree to 
which the charge is confined within narrow stripes. By 
reason of their greater absolute intensities, the magnetic 
harmonics are best suited for quantitative analysis. 

For the e = ^ phase found at 80 K, there are 9 in- 
dependent harmonics. Figure [)] illustrates the spectrum 
of higher order magnetic and charge satellites found in a 
scan along [h, 0, 1) for < h < 1. Because of the weak- 
ness of the higher harmonics, the intensity is plotted on 
a logarithmic scale. The n th order satellites of magnetic 
origin are denoted by the subscript m, while those due to 
charge order are indicated by the subscript c. An oxygen 
ordering peak at (|, 0, 1) (labeled O) and an unexplained 
peak close to n m = 3 (labeled *) are also recorded. (The 
latter very weak peak is probably due to a misaligned 
grain.) Because of insufficient cool-down time for this 
special scan, the first magnetic and charge satellites are 
split in two due to partial coexistence of the and | 
phases. Higher order contributions of the | phase could 
not be observed in this scan. 

In order to allow a quantitative analysis, separate mea- 
surements (following a more careful cool down) were per- 
formed to obtain integrated intensities for all indepen- 
dent magnetic harmonics. These were done at 80 K and 
100 K to characterize the yg and | phases, respectively. 
All peaks in the range Q = (h, 0, 1) with — 1 < h < 1 were 
measured, and values at +h and —h were averaged. The 
results, normalized to the i nten sity of the first harmonic, 



are listed in Tables VI and VII 



To model the harmonic intensities, we must start with 
a specific model of the spin density within a unit cell of 
the stripe pattern, such as the models shown in Sec. IV. B. 
Within such a one-dimensional unit cell, if we number the 
sites with an integer index m and denote the spin on the 
m th site as S m , the discrete Fourier transform of the spin 
arrangement corresponding to the n th harmonic is then 



1 

m=0 



— iixm — i27rnern 



(28) 



where N = 36 (14) for e = y| (§). The intensity of the 
n th -order satellite is given by 
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\F n 



(29) respectively. 



where F n is given by Eq. fl27|) . Because the measure- 
ments span a very small range of Q, the variation of the 
magnetic form factor and the Debye- Waller factor can be 
neglected. 

We first consider the | case, because it provides the op- 
portunity to distinguish between three different models: 
(a) equally-spaced stripes centered alternately on Ni and 
O rows, (b) stripes with alternating spacings centered 
only on Ni rows, or (c) stripes centered only on O rows. 
These three possibilities are illustrated in Eqs. (p~5|)-(17), 



respectively. In each case, all spins are drawn with the 
same magnitude (unity) , suggesting quite narrow charge 
stripes. To allow for some smoothing of the spin mod- 
ulation, we treat the relative magnitude of a spin im- 
mediately adjacent to an O-centered stripe as a fitting 
parameter. Since no such sites occur in model (b), in 
that case alone we take the magnitude of the spins next 
to Ni-centered stripes as the fitting parameter. The in- 
tensities calc ulated with each of the models are presented 
in Table VII , and the fitted spin in each model is listed in 



Table VIII. Model (a), with equally-spaced stripes, gives 



very good agreement with experiment, and the other two 
models are ruled out rather convincingly. As a further 
test of model (a), we allowed a second parameter, the size 
of the spin adjacent to a Ni-centered stripe, to vary, and 
found that it decreased only slightly, from 1.0 to 0.94, 
while the size of the spin next to an O-centered stripe 
remained unchanged. 

To model the results for the jg phase we considered 
only a model similar to (a), with both Ni- and O-centered 
stripes, although in this case they are not all equally 
spaced. The fit gives the same value as in the | case for 
the spin magnitude next to an O-centered stripe (see Ta- 
ble VIII). The calculated harmonic intensities are listed 
in Table VI. The R factor is somewhat worse than in the 



j case, but we also have a substantially larger number 
of harmonics to fit. The quality of the fit seems quite 
reasonable given the relative complexity of the structure 
and the possible systematic errors in measuring the very 
weak harmonics. 

Given that the analysis of the magnetic harmonics in- 
dicates fairly narrow charge stripes, it is of interest to 
study the charge-order harmonics. The information ob- 
tained about the charge modulation is still indirect, be- 
cause the neutrons are only directly sensitive to the in- 
duced atomic displacements. Furthermore, quantitative 
modeling is difficult because of the strong Q dependence 
of the structure factor, even for a single layer Neverthe- 
less, an attempt was made to detect higher harmonics. 
Since the lowest-order superlattice peaks for charge or- 
der correspond to n = 2, the next harmonic is n = 4. 
The strongest n = 2 peak within the accessible Q range 
is at (4 — 2e,0, 1). Searching in the nearby region we 
found finite n — 4 peaks at (4 — 4e, 0, 1) with I = and 
2. At T — 8 K the intensities of the latter 2 peaks, rela- 
tive to the strong n = 2 peak, are 1.4(7)% and 2.2(7)%, 



VII. DISCUSSION 

In the Landau model for coupled charge- and spin- 
density-wave order parameters ,L-3 there are two possible 
ordering scenarios. If, on cooling, the first ordering tran- 
sition is driven by the pure charge term in the free energy, 
then spin ordering will tend to occur at a lower tempera- 
ture; otherwise, the spin and charge order will appear si- 
multaneously. Our present results for the S — 0.133 phase 
of La2Ni04+5 show that charge order appears first, con- 
sistent with the behavior found in La2- a; Sr :r Ni04 (Refjj, 
@, |, H) and in Lai. 6 _ x Ndo.4Sr x Cu04 (Ref. |37|, |$. 
Thus, it is the free energy associated with the charge 
alone that drives the ordering. Of course, one must keep 
in mind that the charge stripes would not be likely to 
form and order-in the absence of the magnetic moments 
of the Ni ionsEj In the charge-ordered phase at T > T m , 
the Ni spins remain strongly correlated, as demonstrated 
by inelastic, neutron scattering measurements reported 
elsewhere. EJ 

The exponential decrease of the charge-order peak in- 
tensity at high temperatures is unusual, but appears to 
be associated with a commensurability with the super- 
lattice formed by the interstitial oxygens. This coupling, 
together with the simultaneous appearance of interstitial 
order and charge-stripe order as a function of 5, leads 
one to wonder whether the energy associated with the 
charge and spin correlations might influence the ordering 
of the interstitial oxygens. Perhaps this connection could 
be tested in a future experiment. 

The first-order transition at T m and the concomi- 
tant jump in / c h are consistent with the generic Landau 
modelB The main feature not explained by that model 
is the temperature dependence of e. The similarity of the 
series of lock-in transitions to what is observed in solu- 
tions o f anis otropic next-nearest-neighbor Ising (ANNNI) 
modelscSEj suggests the importance of competing inter- 
actions. (Of course, it is likely that competing interac- 



tions are also r£3PQ] 
in this systemO'ol 



. e for the very existence of stripes 
3) One relevant factor is the long- 
range part of the Coulomb interaction, which favors a 
uniform, distribution of charge, and hence a large value 
of c.BejO'O A competing factor is the superexchange 
energy associated with antiferromagnetic Ni spins. This 
interaction favors w4de,magnetic domains, and therefore 
a small value of e.E5c3 The competition between these 
energies will depend on the relative magnitudes of the 
order parameters p and S, resulting in a temperature- 
dependent e. The observation^ of a very similar tempera- 
ture dependence (without the lock-ins) in La2- a; Sr a ;Ni04 
with x — 0.225, where the dopant ions are randomly 
distributed, demonstrates that this behavior does not re- 
quire well ordered dopants. 

The lock-in transitions provide evidence for some cou- 







pling to the lattice. We have explained how the values of 
e observed at lock-in plateaus can be understood in terms 
of long-period commensurate structures. Only commen- 
surate structures with relatively short periods can be 
stabilized at finite temperature. Between the plateaus 
entropy will tend to favor incommensurate structures .E3 
The values of epsilon corresponding to observed com- 
mensurate structures may at first appear to be some- 
what arbitrary rational fractions; however, it is inter- 
esting, to. note that they form the first levels of a Farey 
treeJl3'C3 A similar correspondence was ©beerved previ- 
ously by Shimomura, Hamaya, and FujicZI in an X-ray 
diffraction study of commensurate-incommensurate tran- 
sitions in the system [N(CH3)4]2MCi4, where M = Zn, 
Fe, and Mn. A Farey tree is constructed by evaluating a 
sequence of Farey mediants. Given two rational fractions, 
n'/n and m'/m, the Farey mediant is (n' +m')/(n + m). 
Starting with the fractions j and | one obtains the tree 
shown in Fig. [Io| The experimentally observed values 
of e are marked. Farey-tree structures have also been 
noted in studies of the two-dimensional Falicov-Kimball 
model in regions of .-the phase diagram dominated by 
stripe configurations .c3 

It was shown previouslyQ in a neutron diffraction study 
of stripe order in a S — 0.125 sample that the ordered 
Ni moments are quite large, with a maximum amplitude 
that is > 80% of that observed in undoped La2Ni04. 
Here we have shown that the observed intensities of mag- 
netic harmonics are consistent with fairly narrow domain 
walls. These results, together with the poor conductiv- 
ity of the nickelates in general, suggest that the dopant- 
induced holes are strongly localized within relatively nar- 
row charge stripes. The weakness of the higher harmon- 
ics relative to the first harmonic peak occurs because of 
the small number of atomic sites between domain walls: 
the difference in harmonic intensity between a square 
wave and a sine wave is small when there are only a 
few lattice sites per period. It follows that the absence of 
strong higher harmonics in magnetic scattering studies of 
copper-oxide superconductors such as L^-^Sr^CuO^ is 
not inconsistent with the existence of stripe correlations 
in those materials. 

It is interesting to consider the change in character of 
the charge_stripes with temperature. As we have shown 
elsewhere,t3 the charge stripes in the e = | phase, above 
T m , are centered on oxygen rows. From the analysis of 
the harmonic intensities, we have shown here that, in the 
e = | phase, the charge stripes are half O-centered and 
half Ni-centered. Continuing to lower temperature, the 
domain walls become dominantly Ni-centered. In fact, 
this shift of the stripes from O- to Ni-centered is corre- 
lated with the evolution of the magnetic order parameter. 
We have observed that, at least near the transition, the 
intensity of the first magnetic harmonic is proportional to 
the square of the quantity q defined by Eq. (|^). One can 
easily show that 3q is equal to the density of Ni-centered 
stripes. The significance of this curious correlation is not 
clear. 



Along with the shift in lattice alignment, we have ar- 
gued that the charge density within the stripes must vary 
with temperature. It follows that the insulating nature 
of the material cannot be explained simplv— by invoking 
a Peierls transition in ID half-filled stripes.E-3 Instead, it 
appears that the electronic localization must be under- 
stood in terms of interactions in the direction transverse 
to the stripes. While electron-phonon interactions are 
often discussed, the temperature dependence of e makes 
it clear that purely electronic effects are more impor- 
tant than |60|trt)ling to the lattice. It has been suggested 
previouslyc3S that "magnetic confinement" effects are 
relatively strong in nickelates, especially when compared 
to cuprates. Perhaps this effect, associated with the size 
of the transition-metal moments, is the most important 
feature distinguishing the nickelates from the cuprates. 

VIII. SUMMARY 

We have reported a detailed study of stripe ordering 
in a single crystal of La2Ni04+a with S = jt-. On cool- 
ing, charge order occurs before magnetic order, although 
there is a substantial jump in the charge order parameter 
at the first-order magnetic transition. No clear charge- 
order transition has been observed, as the intensity of a 
charge-order superlattice peak decays exponentially with 
temperature. Below T m , the value of e decreases with 
temperature, exhibiting several lock-in transitions. A 
model in which each charge stripe is centered on either a 
row of Ni or O ions has been shown to be consistent with 
the ^-dependence of the magnetic peak intensities and 
with the intensities of higher-order magnetic harmonic 
peaks. Modelling of the latter also indicates relatively 
narrow magnetic domain walls. With decreasing tem- 
perature the charge stripes are initially all centered on O 
rows, and then begin to shift to Ni rows below T m , with 
Ni-centered rows dominating at low temperature. 
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FIG. 1. (a) Diagram of the (hOl) zone, showing positions 
of diffraction peaks, with indexing based on space group 
F4/mmm for the average nuclear structure. Solid circles, fun- 
damental Bragg peaks; solid triangles, magnetic superlattice 
peaks; solid squares, nuclear superlattice peaks correspond- 
ing to charge order; open squares, allowed but unobserved 
charge-order peak positions. Oxygen ordering peaks are ex- 
cluded, (b) Model of spin and hole ordering in a MO2 plane, 
as discusssed in the text; circles indicate positions of holes 
(Ni centered), arrows indicate correlated Ni moments. The 
dashed line outlines a unit cell, and double lines mark posi- 
tions of domain walls. 



11 



La NiO, 8=^ 

2 4+5 ' 15 



1.0 

0.8 h 

m° 0.6 
i— i 

0.4 

0.2 

0.0 
3.45 

3.40 

3.35 h 

3.30 

3.25 



5-H 



10° 



if 

li 



io- 



I 1 I 1 I 1 I 1 I 1 I 




■ ■■■■■ ■ 



120 160 200 
Temperature 



1-8-+- 1 r'*T*T*T»T«t« i 



Q = (h,0,l) 
• h = 4-2e charge 

\ □ h = 3+e magnetic 



□ □□□□□□n 







250 



50 100 150 200 

Temperature (K) 

FIG. 2. Temperature dependences of a magnetic peal 
at (3 + e, 0, 1) (open squares) and a charge-order peak a1 
(4 — 2e, 0, 1) (filled circles). The upper panel shows integratec 
intensities normalized with respect to 10 K, the lower give; 
the peak position in h. Inset: logarithm of normalized inten- 
sity vs. temperature for the charge-order peak above 110 K. 
Line through points is a linear fit, as discussed in the text. 
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FIG. 3. Elastic scans along Q = (h, 0, 1). Scans were mea- 
sured at temperatures of 110 K (a), 110.5 K (b), and 111 K 
(c). 



12 



La.NiO, 8 = — 

2 4+5 ' 15 



La NiO„ s , 8 = 

2 4+8 15 




0.0 1.0 2.0 3.0 

(1/3-e) 2 (xl0~ 3 ) 

FIG. 4. Normalized intensity of a magnetic peak at 
(3 + e,0, 1) vs. (1/3-e) 2 ~ for 80 K< T < T m . Line through 
points is a linear fit, as discussed in the text. 
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FIG. 5. Intensity of the (4 — 2e, 0, 1) charge-order peak, 
plotted in the form AI^ 2 (T) and normalized to the 10-K 
value, vs. the normalized (3 + e, 0, 1) intensity. 
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FIG. 7. Temperature dependences of third harmonic mag- 
netic peaks at (1 — 3e, 0, 1). At most temperatures two com- 
ponents with different e are coexisting. In all panels, the 
(usually, see text) dominant component is represented by the 
filled circles and the secondary by open triangles. The top 
panel shows the integrated intensity normalized to the total 
intensity from both components at 10 K. The middle panel 
and inset show e and the bottom panel the peak width, di- 
vided by 3 (in units of 2-k/o). 
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FIG. 8. Comparison of the structure factors Si (I) for the 
experimental data (corrected for Af 2 ) and least-squares fit 
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are labeled with subscript m, the ones due to charge order 
with c. An oxygen ordering peak is labeled O and one of 
undefined origin with (*). 
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FIG. 10. Farey tree composition rule for wave vectors (den- 
sities) between 1/4 and 1/3. The observed rational values for 
e are marked by the symbols O an d A for La2Ni04.i33 and 
La2Ni04.i25, respectively. 
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TABLE I. Some of the more stable stripe configurations in 
the range \ < e < i assuming a single type of charge stripe. 
The notation is explained in the text. 



n 


m 




e 




Configuration 








fraction 


decimal 




2 


1 




3 

n 


0.2727 


(443) 


5 


3 




8 

29 


0.2759 


(44344343) 


3 


2 




5 

18 


0.2778 


(44343} 


1 


1 




2 
7 


0.2857 


(43> 


2 


3 




5 

17 


0.2941 


(43343} 




TABLE II. 


Some of the more stable stripe configurations 


in 


the range 


\< 


e < | assuming two types of stripes. The 


notation is explained in the text. 




n 


V 


rn 


e 




Configuration 








fraction 


decimal 




1 


1 





3 
11 


0.2727 




2 


3 





8 

29 


0.2759 


(41114111) 

\^2 2 2^2 2 2/ 


1 


2 





5 
18 


0.2778 


(4llll\ 
\^2 2 2 2 / 





1 





2 
7 


0.2857 


(ID 





2 


1 


5 

17 


0.2941 


/ 7 7 7 7o\ 
\ 2 2 2 2 °/ 



TABLE III. Examples of \F\ 2 for both magnetic and charge-order peaks in the case of 
body-centered stacking of stripe-ordered layers. 



e 


|F(l±e,0,0| 2 
I even / odd 


|F(3± 
I even 


e,0,/)| 2 
I odd 


|F(2±2e,0,0| 2 
i even I odd 


|F(4±2e,0,0| 2 
i even I odd 


1 

4 
3 
11 

5 

18 
2 
7 
1 

3 


2 

2 
2 
4 


2 
4 
2 
2 



2 
4 
2 
2 



2 

2 
2 
4 



4 


4 


4 

4 
4 



4 
4 
4 
4 
4 



TABLE IV. Values of the Z-dependent coefficient B n (de- 
scribed in the text) for harmonics n = 1-4. For n = 1, exper- 
imental and fitted values are given. For n = 3, experimental 
values are compared with values calculated using the param- 
eters determined from the n = 1 fit. Only calculated values 
are given for n = 2 and 4. 



T 


Bi 




B 3 




B 2 


Bi 


(K) 


expt. 


fit 


expt. 


calc. 


calc. 


calc. 


80 


-0.13(1) 


-0.132 


0.34(7) 


0.398 


-0.964 


0.841 


100 


-0.15(1) 


-0.153 


0.54(2) 


0.508 


-0.854 


0.508 



1G 



TABLE V. Fitted values ol the occupation fraction a 
for layers with basis vector (ma, 0, §) in the average mag- 
netic unit cell and a, the Ni magnetic form factor constant. 
R w = J2 w n(Iobs - In) 2 / w nllbs\ w n = statistical weight). 



T 


e 


m 


a 


a 


R W 


80 K 


5 

18 


8 


0.123(2) 


3.2(2) 


0.034 


100 K 


2 
7 


3 


0.1182(2) 


3.1(2) 


0.035 



TABLE VI. Relative intensities of the n* order magnetic 
satellite peaks in the Q range (h,0, 1) with —1 < h < 1 for 
e = Both observed and fitted values are listed. 



n 


(/n//l)obs (%) 


(/n//l)calc (%) 


3 


1.70(21) 


1.008 


5 


0.10(10) 


0.168 


7 


0.10(7) 


0.024 


9 


0.001(1) 


0.136 


11 


0.001(1) 


0.014 


13 


0.90(16) 


0.642 


15 


0.10(3) 


0.309 


17 


0.10(5) 


0.121 


Rw 




0.47 


TABLE VII. 


Same as previous table, but for e = |. The 


different models 


are described in the text. 


n (/„//i)obs (%) 


(In/hUlc (%) 




model (a) 


model (b) model (c) 


3 0.81(19) 


0.910 


0.355 1.836 


5 0.04(4) 


0.051 


0.063 0.217 


Rw 


0.13 


0.56 1.42 


TABLE VIII. 


Values for the relative magnitude of the Ni 


spins adjacent to a domain wall obtained in fitting the rela- 



tive intensities of the higher order magnetic satellite peaks. 
Models (a), (b), and (c) are described in the text. 



model 


Fitted reduction factor for (SNi) near DW 




c 7 c 18 


(a) 


0.65(2) 0.66(6) 


(b) 


0.44(4) 


(c) 


0.26(8) 
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